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Abstract 

Charge asymmetry in electron (positron) scattering arises from the interference of the Born 
amphtude and the box- type amphtude corresponding to two virtual photons exchange. It can be 
extracted from electron proton and positron proton scattering experiments, in the same kinematical 
conditions. Considering the virtual photon Compton scattering tensor, which contributes to the 
box-type amplitude, we separate proton and inelastic contributions in the intermediate state and 
parametrize the proton form-factors as the sum of a pure QED term and a strong interaction term. 
Arguments based on analyticity are given in favor of cancellation of contributions from proton 
strong interaction form factors and of inelastic intermediate states in the box type amplitudes. In 
frame of this model, with a realistic expression for nucleon form-factors, numerical estimations are 
given for moderately high energies. 

PACS numbers: 
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I. INTRODUCTION 



The nucleoli structure is traditionally investigated using electromagnetic probes, and as- 
suming that the interaction occurs through the exchange of a virtual photon, which carries 
the momentum transfer, g, from the incident to the scattered lepton. Recently, large at- 
tention was devoted to 27 exchange amplitude both in scattering and annihilation channels 

1 n n fi 

l|, |2|, [SL IJ], in connection with experimental data on electromagnetic proton form factors 
(FFs) [st 

The extraction of box type (two photon exchange amplitude (TPE)) contribution to the 
elastic electron-proton scattering amplitude is one of long standing problems of experimental 
physics. It can be obtained from electron proton and positron proton scattering at the same 
kinematical conditions. A similar information about TPE amplitude in the annihilation 
channel can be obtained from the measurement of the forward-backward asymmetry in 
proton-antiproton production in electron-positron annihilation (and from the time-reversal 
process). 

The theoretical description of TPE amplitude is strongly model dependent. Two reasons 
should be mentioned: the experimental knowledge of nucleon FFs is restricted in a small 
kinematical region, and the contribution of the intermediate hadronic states can be only 
calculated with large uncertainty, the precision of the data being insufficient to constrain 
the models. 

A general approximation for proton electromagnetic form-factors follows the dipole ap- 
proximation: 

GE{q') = = = (1 + QVO.71 GeV2)-^ = -g^ = -t, (1) 

/i 



where /i is the anomalous magnetic moment of proton. However, recent experiments 
showed a deviation of the proton electric FF from this prescription, when measured following 
the recoil polarization method (gI, which is more precise than the traditional Rosenbluth 
separation [3]. Such deviation was tentatively explained, advocating the presence of two 
photon exchange contribution. 

The motivation of this paper is to perform the calculation of charge odd correlation 
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in the process of electron-proton scattering in frame of an analytical model (AM), free 
from uncertainties connected with inelastic hadronic state in intermediate state of the TPE 
amplitude. In frame of this model it is possible to show that there is a compensation between 
the effects due to strong interaction FFs and those due to the inelastic intermediate states, 
within an accuracy discussed below. 

Our paper is organized as follows. In part II, we introduce a new decomposition of proton 
FFs, separating QED and strong interaction contributions. In section III we formulate the 
analytical model and calculate the contribution of the QED part of FFs. In Section IV the 
resulting expression for the asymmetry is obtained. In Section V we present the results of 
numerical integration for asymmetries and in the Conclusions we estimate the accuracy of 
the obtained results. The Appendix contains some details of the calculation. 



II. NEW FORM OF PROTON FORM FACTORS 

The cross-section of elastic ep scattering 

e(pi) + p{p) e{p[) + p{p) (3) 

in Born approximation, in laboratory (Lab) frame {p = (M, 0,0, 0), pi = £"(1,1,0,0)) has 
the form: 

t = q2 ^ 2pip[, s = 2ME, u^--^ -2pp[, 

P P 

Q 

s + i + ii = 0, = 1 + 2(1 + T)tan2-, 

with 

<Tred^TGlt + eGl, Gm^F. + F^, GE^F,-rF2. (5) 

Here 9- is the Lab electron scattering angle. 

In the analysis of the TPE amphtude we consider the electromagnetic interactions in the 
lowest order of perturbation theory. Hadron electromagnetic FFs in the space-like region, 
which parametrize the interaction with the external electromagnetic vertex as 

u{p')[^i{Q'h, + ^qi.r^iq'Mp), q^p'-p,p'^ p' = M^ (6) 



are functions of one kinematical variable, q^. The static value of the Dirac FF of the proton 
(for (5^=0) is unity due to 'QED' origin, i.e. interaction with a point-like quark. Therefore we 
decompose the Dirac FF in two terms, corresponding to different effective internal momenta: 

Fi{q^) = FiQ{q^) + Fu{q^), F^^q^) = F^s^q^), with Fi,(0) = 0, and ^2,(0) = (7) 

where the subscript 's' stands for strong interaction and Fiq corresponds to the 'QED' 
contribution. FiQ{q^) is relevant at very small momenta compared to the typical hadronic 
mass and describes the exchange of gluons with virtuality <^ <^ Q^, where rrig is the 
mass of current quarks. Fiq is a fast decreasing function of Q^, and may be represented by 
an expression of Sudakov type FF: 



FiQiq' 



exp 



4:71 



mi 



(8) 



Contrary to the charged lepton Dirac form-factor, which can be considered constant and 
equal to unity, independently on Q^, the Sudakov parametrization gives a negligible con- 
tribution as deviates from zero. Although the Sudakov parametrization can not, in 



principle, be applied for ^ w?, taking 



ma ~ 



3-5 MeV, the parametrization from ([8]) 



applies starting from the quark mass value. This justifies the following parametrization to 
the hadronic electromagnetic FFs that will be used later on for practical purposes: 



F^q{Q^) = 1,Q^ <m, 



2 2 

~ m- 



F,Q{Q^) = ^,Q^>ml. 



(9) 



In most of the kinematical domain, one can safely neglect the interference between Fiq and 
Fis, as they act in different regions of transferred momenta: F^ = F^q + Ff^. However, 



the derivative at 
interference term: 



IQ 

0, which is connected to the charge radius, r^, is sensitive to the 
d 



Parametrization of space-like form factors in terms of fractional polynomials have already 
been suggested in literature [ol, Q. The present choice of the strong interaction FF is 
consistent with the following parametrization: 



FuiQ' 



/i 







Qf 



F2s{Q' 



1 + Er' h 



(10) 
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where c, d, e,/ can be considered as fitting parameters. The last coefficients of the series 

Cn, dn+3, e„, fn+3 are constrained by the high asymptotic hmit. The low properties: 

Fis(O) = 0, F{ (0) = -rl and -^2^(0) = fj,, where fj, is the magnetic moment, are explicitly 

6 ^ 

taken into account. 

III. FORMULATION OF THE ANALYTICAL MODEL 

Let us discuss now the arguments in favor of mutual cancellation of the terms of order of 

with the contribution of the inelastic hadronic intermediate states, in TPE amplitude. 

The TPE amplitude contains the virtual photon Compton scattering tensor. It can be 
split in two terms, when only strong interaction contributions to Compton amplitude are 
taken into account. One term (the elastic term) is the generalization of the Born term with 
the strong-interaction FFs at the vertexes of the interaction of the virtual photons with 
the hadron. We suppose that the hadron before and after the interaction with the photons 
remains unchanged. The second term (inelastic) corresponds to inelastic channels formed 
by pions and nucleons or the excited states of the nucleon such as the A resonance. 

Denoting the loop momenta in TPE amplitude as k, the loop momenta phase volume can 
be written as: 

d^k = — d'^k\_dseds„, 

where Se = {p~ — k)"^ and Sp = {p + k)"^ are the invariant mass squared of the electron and 
proton blocks, and k±_ is a two-dimensional euclidean vector orthogonal to the momenta of 
the initial particles :k±p- = k±p = 0. Considering both Feynman diagrams, the box and the 
crossed diagram (Fig. [1]), the integration over Se {—oo < Se < oo) can be done calculating 
the residue by in the electron Green function (omitting higher QED corrections). The 
contour for the Sp integration has a Feynman contour form: — oo — iO < Sp < oo + iO (see 
Fig. E]). 

Let now consider the analytic properties of the Compton scattering amplitude with both 
photons ofF-mass shell as a function of the complex variable Sp. The singularities on the 
physical sheet are a pole, located at Sp = M^, which corresponds to a proton in the interme- 
diate state and a series of cuts corresponding to inelastic states of a nucleon accompanied 
by pions, and nucleon-antinucleon pairs. The ffist cut lies at s-j^p = (M -|- m^)^. The left 
cut lies at S2 < — (SMat)^ and corresponds to a ppp state in w-channel. Its contribution, for 
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FIG. 1: Feynman diagrams for two-photon exchange in elastic ep scattering: box diagram (a) and 
crossed box diagram (b). 

the case ~ s ^ is suppressed by powers of (M^/Q^)", compared with contributions 
corresponding to singularities of the right cut, since the operator of higher twist become 
relevant |11|. Neglecting the contribution of the left cut the Sp integration contour can be 
closed to the pole and the Sp > s-^p cuts, insuring their mutual cancellation. Similar con- 
siderations were developed in frame of QED The application of this result to forward 
elastic ep scattering amplitude in the high energy limit allows to derive sum rules for the 
strong interaction contribution to the proton FF. The contribution of the left cut, in this 
case, has to be precisely calculated. This will be published elsewhere. 

Based on this hypothesis (which was proved rigorously in frame of QED and here taken 
as an assumption) we consider only the QED part of the amplitude corresponding to one 
nucleon (non-excited) state in the hadronic block and omit the pure strong interaction 
contributions. 

Proton FFs enter in the box amplitude in a form which can be schematically written as 



m2(e)(p) {k) (k) ' ^^^^ 

(A:) = P-A^ Ck) = ik-qr-\^ {e) = {k - p,f - ml {p) = {k + pf - M\ 

where we extract the 'QED' part and do not distinguish Dirac and Pauli form factors. This 
integrand expression can be rearranged as 

F,Q{k^)F,Q{¥) F,{k')F,m , F,Q{e)F,{k') , F,{k')F,Q{k') ...^ 

{km mk) {km {km ' ^ ^ 




(2NN) 



(a) 




FIG. 2: Illustration of singularities along the S2 real axis with the open contour C (a), and with 
the contour C closed (b). LC stays for large circle contribution. 



The first term cancels due to the assumption about the 'QED' contribution to the Dirac 
form-factor iQ. The second one can be omitted due to the hypothesis of cancellation in 
strong interaction. The last two terms coincide and can be written as: 



-9{ml-\k'\). 



(13) 



i7r^(e)(p) (k) 

The relevant integral is calculated in the Appendix. 

The virtual photon emission contribution to the cross section has a form (according to 

Eq. my- 



a 



day 

dO 27rt2M2p2' 



(14) 



with: 



d^fc 1 St 
m'^ {k) {p) 



9iml-\k'\), 



(15) 
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and 



(e) = (fc + Pi) — m 



Sj' 



^Tr{p + M)r,{~q){p' + M)T,{k){p + k + M)T ^{q - k), 



(16) 



Using the formulas given in Appendix, the virtual photon contribution to the differential 
cross section with two photon exchange, da^, can be written as: 



, „ , , 2a , / 2EM 1 , , , 
da^ = — In p I In 2 ' daBom- 



(17) 



The IR divergence from virtual photon emission contribution is, as usually, canceled when 
summing the contribution from emission of soft real photons: 



dam doBbox] .odd ^ da^r , da'°^^^^ 
dn dn \ dfi dfi ' 



(18) 



The quantity 5°^^^ was considered in 



an 



S^odd 
"soft 



_2 / d^^Pi_ _ pi p' 



p 



2a 

TC 



IQtt^ J io p[k pik p'k pk' 

- 1 2pAE 
m — m 



5o,a;<AE 



A 



+ Ina;lnp + Li2 ( 1 - — ) - Lig (l - - 

\ px J \ X 



X 



with A, and AE respectively the mass and the maximal energy of the soft photon. 



(19) 



IV. RESULTS AND DISCUSSION 



The final result for the asymmetry A°'^'^ 

j^odd 



is: 



2a 

TT 



1 



P 



' 1, (2AE)2 5, 2 , , 
In - In — — — In p + In a; In p + L12 1 

p ME 2 \ px 

s 



Li2 (1 



(20) 
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The finite part of tlie asymmetry is calculated for different values of and Q. The results 
are plotted in Fig. [3l as a function of Q. The asymmetry vanishes for = and reaches 
the largest values for 9 = -k. It is measurable, of the order of 5 %. The largest contribution 
to the asymmetry is given by the first term in Eq. fl2U]) . which depends on the soft photon 
energy. The calculation is done for /S.E = O.OIE. 

In Ref. [l^, the approach used to calculate TPE has to be considered as a model: one of 
the exchanged photon is quasi real. The contribution of inelastic intermediate states is also 
ignored. The comparison with the results from Ref. [l^, is shown in Fig. HJ where 

the absolute difference 

AA = A^^^ - A""^ = ^ (^In 1 In ^ + 1 In^ (21) 

between the two calculation is shown. The present results are in general smaller, except at 
small angles, where they are comparable. The main difference is due to the term related to 
ln(Q^/s) which gives a different AE dependent contribution. 

V. CONCLUSIONS 



Charge asymmetry in electron proton elastic scattering contains essential information 
on the contribution of the real part of the 27 exchange to the reaction amplitude. This 
amplitude can shed light on Compton scattering of virtual photons on proton. It contains 
a part corresponding to proton intermediate state, which carries the information on proton 
FFs. Another term corresponds to excited nucleon states and inelastic states such as Nit, 
N2'K, NNN. Their theoretical investigation is strongly model dependent. 

Similar effects of charge and angular asymmetries can also be due to Z-boson exchange, 
but such contribution is small for moderate-high energy colliders. The ratio of corresponding 
contributions can be evaluated as: ~ {TrgvgAs)/{aMz) < 5 ■ 10^'^ for s < 10 GeV^ {gv {gA} 
is the vector (axial) coupling constant of the Z boson to fermion). 

The analytical calculation of 27 amplitude with FFs encounters mathematical difficulties. 
In Ref. ^5|] the results for the box amplitude with arbitrary FFs was investigated. Other 
works use different approaches to include FFs, and the results are quantitatively different 
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FIG. 3: Charge odd correlation in electron-positron scattering as a function of the scattering angle 
e, for Q2 =1 GeV2 (thin sohd line), 3 GeV^ (thick soHd hne), 5 GeV^ (dashed line), 7 GeV^ (dotted 
line), 9 GeV^ (dash-dotted line). The calculation corresponds to AE'=0.01 E (see Eq. [20]) . 



The analytical model presented in this paper is based on two main assumptions: - the 
separation of the nucleon FFs in two terms, one of which corresponds to the QED contribu- 
tion, which contributes explicitly to a narrow region of transferred momenta, close to zero. 
- the mutual compensation of the strong interaction contributions to the TPE amplitude, 
arising from nucleon form factors and from inelastic states. The last assumption, which has 
been proved in QED and which holds for zero scattering angle amplitude, has to be consid- 
ered as an approximation when applied to the high energy limit of large angle scattering. 
The tendency of cancellation of elastic proton state and the A resonance in the real part of 
the TPE amplitude was previously noted in literature (0, Q). 

The numerical results obtained here are in agreement with our previous calculation 

m 

and confirm our previous conclusion that two photon contribution can not be responsible for 
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FIG. 4: Difference between tlie present asymmetry and the calculation [14], A^, Eq. [2T1 Notations 
as in Fig. [3l 



the discrepancy in recent FFs measurements. A more plausible explanation is that higher 
order radiative corrections should be taken into account in the leptonic vertex j4|. 

Other works [16|] devote attention to the excited intermediate states as A and A^* reso- 
nances, introducing additional uncertainties. In our approach, excited states should not be 
included, as they correspond to poles in the second physical sheet. Only contribute nirN 
states, with any number of pions. 

Our main assumption about the compensation of pure strong interaction induced contri- 
butions to FFs and inelastic channels allows us to avoid additional uncertainty connected 
with inelastic channels. Experiments measuring charge-odd observables in ep scattering will 
be critical for the verification of the validity of our model. 

The numerical results show that charge-odd correlations are of the order of few percent, 
in the kinematical region considered here. Such value is expected to be larger at larger 
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values and could be measured in very precise experiments, at present facilities. 
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VII. APPENDIX 

The calculation of the integral 

SMb"'"" - I'^l' '''' 

is performed using the Sudakov's momentum parametrization. For this aim we build two 
(almost) light-like vectors 

Pi = Pi-p — ,P = P-Pi ,s = 2ppi,Pi = 0;p^ = 0, 

s s 

2pp = M2;2pipi = m^. (23) 
The integration is done over the momentum variables: 

k — ap + Ppi + k±, k±p — k±pi — 0,k^ — sa(5 — k'^, 
and the phase volume is : 

df'k — -dad(5d^k. 

The analysis of the location of the poles of the denominators in a and /3 planes: 

{k) = A;^ - = sa(3 - P - + iO; 

(e) = k — pi)^ — n? -|- iO = saji — }? — sa — m^P + iO; 

(p) = (^k + py-M^ + iO = sai3-P + s(3 + M^a + iQ, (24) 

leads to two regions of non zero contribution < a, (3 < 1,-1 < a, (3 < 0, with equal 

contributions. We obtain 
11^ 

r r fiu^ 1 1 1 
1 = 2 da df3-. ^ — ^ 9iml - (25) 
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the real part of / (which is relevant for us) can be extracted using the identity l/(x + iO) 
V{l/x) — in5{x), and performing the integration on k"^: 



dz5{z — a) = 9(a), a = sajS — A^, > a > 0. 



The final answer is 

Rel = 2 



- In^ h In In — In^ 

2 mM mM 2 m 2 



It is interesting to compare this result with the same integral without cut on |A;^ 
d^'k 1 



J = Re 



In 



^ s . M , s , mM 47r2 

2 - In^ — + 2 In In 



(26) 



(27) 



ivr^ {k){e){p) mM m mM 3 

Using similar expressions for another integral of this type, it is straightforward to recover 
Eq. (1201). 
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